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Abstract 

We consider first-passage percolation on a ladder, i.e. the graph 
N x {0, 1} where nodes at distance 1 are joined by an edge, and the 
times are exponentially i.i.d. with mean 1. We find an appropriate 
Markov chain to calculate an explicit expression for the time constant 
whose numerical value is sa 0.6827. This time constant is the long-term 
average inverse speed of the process. We also calculate the average 
residual time. 



1 Introduction 

Consider a graph G with vertex (node) set V and edge set E C V x V. An 
(undirected) edge e = (v,v') = (v',v) joins vertex v and v' . A path ir(v,v') 
between v and v' , if it exists, is an alternating sequence of vertexes and 
edges (vq, e\, v i, . . . , e n , v n ) such that e$ = (vi-\,Vi) for i = 1, . . . ,n, urj = v 
and v n = v' . 

Associate with each edge e a non-negative random variable £ e . Let the 
time of a path tt(v, v') be 

Ttt(v, v') = \ £ e and let T(v, v') = inf Tir(v, v') 

r-i /^ 7r(«,ti') 
e£ir(v,v' ) 

be the shortest time of any path between v and v'. This may be called the 
first passage time of v' (with respect to v) and is the subject of investigation 
in first passage percolation, see e.g. |S W78j . 

We will think of this as a model for a contagious disease and refer to the 
first passage time of a node as the time of infection. Usually, we think of 
some node v (or more generally a set of nodes) as infected at time zero. The 
times T(w) = T(v,w), w G V, lets us know when nodes w G V are infected. 

One natural question is; how fast is the spread of the infection? The 
present paper considers the equivalent question of the inverse speed of per- 
colation on a ladder, i.e. how long does it take for the infection to spread one 
more step up the ladder, averaged over time. We also calculate the average 
inverse speed in another sense; namely the average residual time. 
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2 First-passage percolation on the ladder 



First-passage percolation on a ladder has previously been studied in |FG806| . 
which gives a method of calculating the time constant when the times as- 
sociated with edges have a discrete distribution, as well as a method for 
getting arbitrarily good bounds for the same quantity when the distribution 
is continuous (and well behaved). Our paper gives an exact formula for the 
time constant only when the times are exponential and the method we use 
relies heavily on the special properties of this distribution. 

It has been brought to our attention that during the referee period of 
the present paper |Sch09j has emerged, which also contains the percolation 
rate (inversion of the time constant), although it is achieved through other 
methods. 

Let L be a ladder, by which we mean a graph with vertex set V = 
N x {0, 1} and where the edge set E consists of (v, v') such that v and v' 
are at distance one from each other, see Figure [TJ By N we mean the set of 
non- negative integers {0,1, 2,...}. 

l o o o o o o o 



o o o- -o o o o 6- 

1 2 3 4 5 6 



Figure 1: Part of the ladder L. 

We will consider the model where the nodes (0, 0) and (0, 1) are infected 
at time zero and each edge is associated with a random variable that has 
an exponential distribution with mean 1. Random variables associated to 
different edges are independent. 

For a node v = (x, y) we will say that v is at height x G N and level 
y G {0, 1}. Recall that T[(x,y)] is the infection time of node (x,y). Let 

= sup{x : T[(x, i)]<t}, i = 0, 1, 

i.e. let Nf: and N^ denote the height of the infection at level and 1, 
respectively, at time t. 

We will let Nt = max{iV t , N^} denote the height of the infection, i.e. 
the largest height of any infected node, at time t. 

We aim to calculate the time constant r, defined as the limiting inverse 
percolation rate, 1/r = lim.tNt/t. We will not do this directly but calculate 
the percolation rate at "time infinity", or rather when a process related to 
Nt has reached stationarity. 



2 



More explicitly, we will employ an observation by Aim and Janson |AJj 
that arose in connection with their joint work on one-dimensional lattices 
[AJ90| . namely; it might be fruitful to consider the front Ft = — 
at time t > 0. Front thus means the absolute value of the difference in height 
between the highest infected nodes on each level, see Figure [2 

Now, the process {Ft,t > 0} behaves like a continuous time Markov chain 
on N, with Fq = 0. Suppose that this Markov chain tends to a stationary 
distribution II = (^0,^1,^2, ■ ■ ■) on N. Then at a late time t the process F 
will be in state with probability ttq. From this state there are two possible 
nodes that may be the next infected ones that result in an increase of the 
./V-process, i.e. the intensity is 2 towards a state that increases the height of 
the infection, since these two nodes are associated with two distinct edges 
(and thus two distinct exponentially distributed mean 1 random variables). 
With probability 1 — ttq the F process is in some other state (F > 1) from 
where there is intensity 1 towards a state that increases the height of the 
infection. Thus, knowing the stationary distribution (or rather knowing ttq) 
gives us the percolation rate at a late time as 

2VT + 1(1 -7T ) = 1 + 7T . (1) 

Below we calculate this stationary distribution II. To do so, it turns 
out that we may express any ir n in terms of ttq in the form ir n = a n iiQ — b n 
where both sequences of (positive) coefficients, a n and b n , satisfy a certain 
recursion (Claim 1) and can be expressed in terms of the Bessel functions of 
the first and second kind. We stumbled upon the solution, and may not have 
solved it otherwise, when it was noted that the sequence b n transformed into 
B n = (b n — b n -{)/n is part of sequence A058797, as listed in [SlOTj . which 
is related to the Bessel functions. This relation to Bessel functions was a 
surprise and we see no reason to expect it in this context. This method 
can be generalized to percolation on a more general "ladder-like" graph, see 
|PR| . The sequence b n is also a subsequence of A056921. 
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Figure 2: Infected nodes at time t marked as black. Here Nt = 6 and Ft = 2. 

Now, we start by looking at the intensity matrix Q of the process Ft. 
Consider the case in Figure [5J where Ft = 2. Now, there are two edges 
leading to the node (5, 1) and if this is the next infected, the process will 
end up in state F = 1, i.e. the intensity to this state is 2. There is one edge 
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each to node (6, 1) and (7, 0) and if these are the next ones infected, the 
process will end up in state F = or F = 3, respectively. This argument 
gives us the third row of the intensity matrix Q. We omit the rest of the 
details, but it is not hard to see that 



/ 



Q 



-2 
2 
1 
1 
1 
1 
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1 
1 






1 

-6 
2 



V 



(2) 



{Ft, t > 0} is irreducible, so if we can find a distribution II = (ttq, tt\, 7T2, . . .) 
on N such that HQ = then II is the unique stationary distribution for 
{Ft,t > 0}, by Markov chain theory. For calculations we henceforth assume 
that such a distribution II exists. 

Using the equations arising from columns one to three of Q, the condition 
HQ = and ^ = 1 we get 



tv i = 3ttq — I, 7T2 = ll7ro — 5, and 7r3 = 567Tq — 26. 



(3) 



Also, taking differences of equations arising from adjacent columns of Q, 
from column three and upwards, gives the equations 



= 7T n _ 3 - (n+ l)7T n _ 2 + (n + 3)7T 



n-1 



7T„ 



n > 4. 



(4) 



We set Ti n = CL n 7To — b n , where 01,02,03,61,62,^3 are defined by (j3J) for 
n = 1,2,3 and, for n > 4 by with a n or b n in place of 7r n . Table [T] gives 
the first nine numbers in the sequences {a n } and {b n }. 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


a n 


3 


11 


56 


340 


2395 


19231 


173490 


1737706 


19136803 


b n 


1 


5 


26 


158 


1113 


8937 


80624 


807544 


8893225 



Table 1: The beginning of the sequences {a n } and {6 n }. 
Claim 1. The sequences {a n } and {b n } both satisfy the recursion 

Cn+l C-n C n C n — 1 . . 

c n = — : , n > 2. (5) 

n+1 n 

Proof. We will use induction. First, check initial conditions, 02 = 11 and 
(56 - ll)/3 - (11 - 3)/2 = 11, 6 2 = 5 and (26 - 5)/3 - (5 - l)/2 = 5, so 
equation dSJ holds for n = 2. Rearranging (JSJ gives 

c n +i - c n = in + 1) ( C n + -[Cn - C n -l] J , n > 2, (6) 
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and this relation is what we will prove. So assume Q holds for some n > 2. 
We aim to show that it also holds for n + 1. So, as n + 2 > 4, we can use 
(H]). Adding ir n — Tt n -\ to both sides of, and substituting n + 2 for n in (jl]) 
gives 

Cn+2 - c n+ i = c n _i - (n + 3)c n + (n + 5)c n+i - c n+ i 

= (n + 2) ^c n+ i + — [2c n+ i - (n + 2)c n - (c w - c^_i) ] ^ . (7) 

~i*r 

From ©, and the induction hypothesis, we get c n — c n _i = (nc n+ i — n(n + 
2)c n )/(n + 1), which we apply on (*) above, and 

17} = (n + 2) (c n+1 + -ij [( 2 - ^) c n+1 - (n + 2 - ^) c n 



□ 



(n + 2) c n+1 H — [c„+i - c n 

n + 1 



Now set C n = (c n — c n _i)/n. Then, if c n satisfies he. c n = C n+ i — C n , 
this implies both 

C n +i = c n + C n (8) 
and C n = C n+ \ — c n . Substitute n — 1 for n in the last equation to get 



Cn-i — C n — C n -\. 



(9) 



Then adding (JSJ) and @ and replacing c n — c n _i with nC n gives us the 
following. 

Claim 2. A n = (a n — a n _i)/n and B n = (b n — & n _i)/n froi/i satisfy the 
recursion 

C n+ i+C n _i = (n + 2)C n , n>3. (10) 



2.1 Bessel functions 

The functions 



and 



7T 



n-1 



2( 7 + ln-J Jn(^)-^ 



(n — k — 1)! fx\ 2k ~ n 



k=0 



k\ 



\) 



2 



fc=0 



k\{n + k)\ \2 



5 



where n 6 N, q m = X^mi J ari( ^ 7 = 0.577 ... is Euler's constant, are known 
as Bessel functions of first and second kind. They both satisfy the recursion 

C n+l {x) + C n _i(s) = —C n {x). (11) 

x 

Define the function 

T(n,m) = vr[J n (2)y m (2) - J m (2)Y n (2)]. (12) 

This function inherits recursion (jlip in parameter n with x = 2, because for 
m fixed T is a linear combination of J n (2) and Y n (2), so that 

Y(n + 1, m) + Y(n — 1, m) = nT(n, m). (13) 

Claim 3. For a// integers n and m the function T(n, m) is integer valued. 

Proof. Clearly T(m, m) = so it suffices to show that T(m+ 1, m) is integer 
valued, since recursion (fT3j) takes care of all other values. 

Two important relations for the Bessel functions are needed. The first 
is that 

xJ' n (x) = nJ n (x) — xJ n+ i(x), Vn. (14) 

and the same holds with Y n replacing J n , see 3.2(4), p. 45 and 3.56(4), p. 
66 of [Wat22j, resp. The second is a fact relating to the Wronskian, namely 
that 

— = J n {x)Y^x) - j' n {x)Y n (x), Vn, x ^ 0, (15) 

7TX 

see 3.63(1), p. 76 of |Wat22] . By dividing both sides of equation (fT4"j) by 
x = 2 and using the resulting expression for J' n {2) and Y^(2), resp., in (fT5|) 
results in 



- = Jn(2) 
7T 



2 



n 



-Y n {2) - Y n+1 (2) - -J n (2)- J„+i(2) y n (2) 



2 

J n+ i(2)y n (2) - J n (2)y n+1 (2) = ir(n + 1, n), 

7T 



the last equality by the definition (|12|) of T. Hence, T(n + l,n) = 1. □ 

Now we can evaluate the function Y(n, m) for different values of n and m 
to find that B n = T(n + 2,0), n = 2,3, and A n = 2T(n + 2, 3) + T(n + 2, 0), 
n = 2, 3, see Table [2 Since T(n + 2, m) satisfies recursion (|10p in n this is 
enough to know that they agree for all n > 2. 

By relation 0, definition (fT2|) and the above expression for A n and B n 
we get 

6 n = T(n + 3,0) - T(?i + 2,0) and 

a n = 2[T(n + 3, 3) - T(n + 2, 3)] + T(n + 3, 0) - T(n + 2, 0). (16) 
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n 


1 


2 


3 


4 


5 


6 


7 


b n 


1 


5 


26 


158 


1113 


8937 


80624 


B n = {b n ~ b n -i)/n 




2 


7 


33 


191 


1304 


10241 


T(n,0) 


1 


1 


1 


2 


7 


33 


191 


a n 


3 


11 


56 


340 


2395 


19231 


173490 


A n = {a n - a„_i)/n 




4 


15 


71 


411 


2806 


22037 


2T(n,3) + T(n,0) 


-3 


-1 


1 


4 


15 


71 


411 



Table 2: Connection between parameters a n ,b n and T(n,m). 
Henceforth we abbreviate J n = J n {2) and Y n = Y n {2), hence e.g. 

We need asymptotics for a n and 6 n . 
Claim 4. 

lim ? ^om = J o lim Q " = 2J 3 + J . (17) 

Proof. |XS64] lists the following asymptotic relations, for fixed x and n 
tending to infinity, 

Jn{ x ) 7= (?) and K„(x) W — ( — ) , 

hence, by the well-known Stirling formula n! ~ \f2~Kn n+1 / 2 e~ n we conclude 
that J n ~ 1/n! and Y n ~ — (n — l)!/vr. Thus 

1 (n - 1)! 

T(n,m) ~7r— y m + 7rJ m , J m • (n - 1)!, 

n! 7r 

from which (|17p follows. □ 

Now we are ready to find an expression for ttq. Since ir n — > as n — > oo 
we get 

no = lim ZEll^L = J ° = 0.4647184275 . . . (18) 

n->oo a n 2 J3 + Jq 
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This in turn gives us 
7T n = a n 7r - b n 

= 2[T(n + 3, 3) - T(n + 2, 3)]tt + [T(n + 3, 0) - T(n + 2, 0)][tt - 1] 
2vr 



2 J 3 + Jo 



[Jn+3^3 — J3^n+3 _ Jn+2^3 + J 3 ^n+2] J( 







— [Jn+3^0 — Jo^n+3 _ Jn+2^0 + Jo^n+2] J3 

2ir(J 3 Y -J Y 3 ) 2T(3,0) 

^— (J n+2 - J n+3 ) - ^^(Jn+2 " Jn +3 ) 

(J n+2 - J n+3 ), n > 1. (19) 



2 J 3 + Jo 



Then ^ 7Tj = 1 — 2Jat +3 / '(2J 3 + Jo) and as J n — )• 0, as n — )• oo, we can 
verify that II is indeed a distribution. 

2.2 The time constant 

Following the discussion before and after ([1]) we get the time constant 

r = — - — = 2 f 3 + J ° = 0.6827250759 . . . (20) 
l + ^o 2J 3 + 2J v ; 

This can be compared with first-passage percolation on N, i.e. the infinite 
path graph (with exponential r.v.'s with mean 1), where the time constant 
necessarily is 1. We may also compare it with 1/2 which is the time constant 
if the rungs of the ladder where associated with random variables all iden- 
tically 0, i.e. there is immediate infection between nodes on different levels 
at distance 1 from each other. 

Note that we started the process by infecting both (0,0) and (0,1). 
Consider the infection of the single node (0, 0) at time 0. Then, the front 
process F will not be a continuous time Markov chain (as before), as the 
behaviour of the front is very different before and after the random time 
Z until there are infected nodes on both levels. Since we are interested in 
stationarity we may simply disregard what happens before Z and view this 
as starting the front process F at time Z, Ft = Fz+t, t > 0. Then F is a 
continuous time Markov chain with intensity matrix given by (J2|), as before, 
albeit started in a random state instead of Fq = 0. Z is dominated by a 
mean 1 exponential random variable (the edge between (0, 0) and (0, 1)). As 
F will converge to the stationary distribution II, so will F. 

Hence, infecting just a single node initially makes no difference to the 
time constant. (In fact, starting with an initial infection of any finite number 
of nodes will not change the time constant.) 

This time constant also gives an upper bound for the time constant 
for first-passage percolation with mean 1 exponential r.v.'s on Z 2 , but not 
stronger than already existing bounds, see |AP02j . 
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2.3 Residual times 



The fact that we can write down the stationary distribution II for the process 
Ft also gives us the opportunity to calculate another statistic for first passage 
percolation on the ladder, namely the average residual time. As far as we 
know, this quantity has not previously been studied for percolation. Define 
the residual time at t as Rt = inf{s : Nt+ S = Nt + 1}, i.e. the time it takes, 
after t, for the infection to spread one more step up the ladder. We are 
interested in T = Ei?t for a late time t. 

One might a priori think that T and r are the same thing. If we were 
looking at percolation on N, which is the same as a Poisson process, then 
they are. This is the waiting time paradox; on one hand: looking at a fixed 
time point increases the probability of choosing a long time interval, but on 
the other: typically half that interval has passed. For a Poisson process, 
these effects balance out. For percolation on a ladder, as shown in (|20|) and 
(|23p below, they do not. 

Now, to calculate T we will assume that F has its stationary distribution, 
so that we may disregard time and put R = R t . Let j n = K[R\F = n]. 
Hence, j n is the expected time it takes for the infection to spread another 
step given that the front is in state n. 

If F = then the residual time is the minimum of two mean 1 exponential 
random variables, so that 70 = 1/2. For n > 1, we write down a recursive 
formula for j n . If F = n then there are n + 2 possible edges, equally 
likely, that the infection might spread along. Only one of these results in an 
increase in the height of the infection, namely the one that gives F = n + 1. 
Two edges result in F = n — 1 and the remaining ones to one of the states 
{0, 1, . . . , n — 2}, respectively. This gives the formula 




(21) 



By writing down the formula for 7 n _i one sees that 



n-2 








which inserted into (|2ip yields the incremental relation 



(n + 2)(7„ - 7„_i) = 7„_i - 7n_ 2 , 



which together with the fact that 71 = |(1 + 2^) = | 
formulas 



in turn gives us the 



7n - ln-l = 



1 



n > 1 



and 




3=0 



(22) 
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Then, as Ft is in state j with probability ttj where it remains on average 
7j amount of time before Nt increases, we get the average residual time T 

as 



T = £ *n7n = 2J — j ( Jo/2 + 2 jT( J n+2 - J n +3)7n ) 
n=0 3 ^ V n=l / 

~+Tq (\ J ° + 2J371 + 2 Yl^n+l ~ 7n ) Jn + 3 ^ 

1 / 1 /I 00 7 \ 

| -J + -J3 + 2V - " +3 N , I = 0.5953444665 . . . (23) 

+ J 1^2 3 ^(n + Sy.) { > 



2J 3 
2Ji 
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